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Variational Mote Carlo Study of Flat Band Ferromagnetism - Application to CeRh3B2 
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A new mechanism for ferromagnetism in CeRh3B2 is proposed on the basis of variational Monte 
Carlo results. In a one-dimensional Anderson lattice where each 4/ electron hybridizes with a ligand 
orbital between neighboring Ce sites, ferromagnetism is stabilized due to a nearly flat band which 
is a mixture of conduction and 4/ electron states. Because of the strong spin-orbit interaction in 4/ 
electron states, and of considerable amount of hybridization in the nearly flat band, the magnetic 
moments from 4/ and conduction electrons tend to cancel each other. The resultant ferromagnetic 
moment becomes smaller as compared with the local 4/ moment, and the Fermi surface in the fer- 
romagnetic ground state is hardly aff'ected by the presence of 4/ electrons. These theoretical results 
are consistent with experimental observations in CeRhsBa by neutron scattering and dHvA effects. 
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1. Introduction 

CeRh3B2 has the highest Curie temperature Tc = 120 K among known Ce compounds. If 
one uses the de Gennes scaling from the transition temperature of 90 K in GdRh3B2,^ the transition 
temperature is estimated as small as IK. Thus the ordinary RKKY interaction cannot be responsible 
for the ferromagnetism in CeRh3B2. The magnetic moment per formula unit of CeRh3B2 is 0.45//b>^ 
which lies in the c plane. It is much reduced from the free Ce^^ ion value 2.14 //b- Because of these 
anomalous features the ferromagnetism in CeRh3B2 has been attracting much interest. 

CeRh3B2 has the hexagonal CeCo3B2 type structure, whose space group is P6/mmm(Dg^). The 
crystal structure of CeRh3B2 is characterized by the lattice parameters a = 5 All A and c - 3.091 
A. The lattice constant along the c axis is remarkably short; shorter than 3.41 A in a-Ce with valence 
close to 4+. Hence a large hybridization between the hgand and 4/ orbitals, which is called the c- 
f hybridization, and a quasi-one dimensional band structure are expected. Previous studies actually 
suggest the quasi-one-dimensional feature along the c-axis.^"^ 

Several models have been proposed to explain the anomalous ferromagnetism. However, mecha- 
nism of ferromagnetism in CeRh3B2 has not yet been identified. For example, an itinerant model of 
ferromagnetism due to Rh 4d can explain small magnetic moment. It was considered as support of 
the model that no magnetic order is observed in CeRu3B2.' However, the model cannot explain the 
fact that Rh Ad bands of LaRh3B2 and CeRh3B2 look essentially the same according to valence-band 
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photoemission,^ while LaRh3B2 shows Pauli paramagnetism and superconducting transition at about 

On the other hand, a 4/ itinerant model can also explain the small ferromagnetic moment.^ It is a 
likely scenario, since the extreme proximity of Ce ions along the c axis may cause some delocalization 
of 4/ charge. Batista et al. have in fact proposed that the ferromagnetism of CeRh3B2 can be explained 
by the periodic Anderson model (PAM).'^ In their work, the high Curie temperature and the small 
magnetic moment are obtained. Recently, de Haas-van Alphen (dHvA) effects have been observed for 
LaRh3B2 and CeRh3B2.^ The dHvA results suggest that 4/ electron is localized, which is in conflict 
with the 4/ itinerant ferromagnetism models. As a result, a localized electron model has also been 
considered as a likely model. However, the small magnetic moment is difficult to be explained by this 
model. 

In present work, we propose a new model for the ferromagnetism in CeRh3B2 with special atten- 
tion to the role of c-f hybridization and the orbital moment of 4/ electrons. According to the band 
calculation,^ the relevant conduction band is regarded as quasi one-dimensional, and is mainly formed 
by molecular orbitals of Rh Ad states lying between Ce sites along the c-axis. We therefore take the 
one-dimensional Anderson lattice with hybridization between the nearest 4/-ligand orbitals. To ana- 
lyze the model, we use the Optimization Variational Monte Carlo(O-VMC) method. We determine the 
ground state phase diagram and analyze the mechanism of ferromagnetism in terms of the effective 
one electron band which is obtained by 0-VMC. On the basis of these results, we discuss the mech- 
anism of the ferromagnetism of CeRh3B2. In this paper we consider properties of the ground state 
only. 

2. Model 

The spin-orbit interaction leads to the total angular momentum 7 = 5/2 for the lowest 4/^ state. 
The six fold degeneracy is split into three Kramers doublets by the crystal field. These states are 
represented in terms of the basis l/^) of 7 = 5/2 as follows: 

|±5/2) = ±^|±3,+l/2)+ ^|±2,±l/2), (2.1) 
|±3/2) = ±^|±2,+l/2)+ y?|±l,±l/2), (2.2) 

1+1/2) = ±^|±I,+I/2)+ ^|0,±l/2), (2.3) 

where \L^, S ^) stands for a basis of L = 3 and 5 = 1/2. 

The crystal-field ground state is taken as |±l/2).^ We keep only the doublet 1+1/2) for explicit 
calculation including hybridization. This simplification is motivated by the high excitation energy of 
220K between |±l/2) and |±3/2),^ which is about two times the transition temperature. When the 4/^ 
states are restricted to |±l/2), the spin-orbit interaction makes the xy plane an easy plane, and the z 
axis the hard axis. To clarify this point, we consider the wave functions where the moment is 
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Fig. 1 . Illustxation of the model in which each / orbital hybridizes with nearest Rh 4d molecular orbitals. 



along the x axis. We obtain 

=^(|l/2> + |-l/2». (2.4) 
V2 

The expectation values of moments are obtained as 

<'P.vtI-/.vI'J'.vT> = 3/2, (l/2|7^|l/2> = 1/2. (2.5) 
Hence the easy axis lies in the xy plane. 

According to dHvA experiment^ and the band structure calculation,^ the Fermi surface of 
LaRh3B2 has the quasi-one-dimensional feature and consists of Rh 4d orbitals. It has been pointed 
out that the origin of the quasi-one-dimensional feature is the strong hybridization between Rh 4d 
bands along the c axis.^'^ We adopt the one dimensional model where the lowest CEF orbital of 4/ 
states hybridizes with a molecular orbital formed by six Rh sites surrounding the c axis threading 
Ce ions. The center of the molecular orbital sits at the mid point between the Ce sites. We take into 
account the nearest neighbor hybridization between Rh 4d and Ce 4/ states. Figure 1 illustrates our 
model. 

For each site we introduce a creation operator d^. for the Rh 4d molecular orbital and /J for the 
Ce CEF states, where cr specifies the the eigenstates of the spin operator S x- We note that (J) 
is anti-parallel to (5> in the Hund-rule ground state. Hence the spin points to the opposite direction of 
the total moment Jx of the 4/ electron. Our model is written as follows: 

-K = ^ tdldj^ + 2] €ddldi^ + 2 efflfi^ 

(ij)cr icr icr 

[ij]cr i 

where (ij) denotes a nearest neighbor Rh orbital pair with hopping energy t, and TjUJ] is the summation 
over nearest 4/ and d sites with hybridization matrix element V. We note that the distance between 
the nearest / and d sites is half of the lattice spacing. The energy of the the 4/ orbital is written as e/, 
and that of the d orbital as ed, which is taken to be 0. Namely we take the origin of the energy at the 
center of the conduction band. We take t = I as the unit of energy. The positive value puts the band 
bottom at the edge of the Brillouin zone, which is consistent with the band structure of LaRh3B2. 

This model has a characteristic band structure in the noninteracting limit U = 0. Namely the lower 
hybridized band becomes completely flat with e/ = V^/t- 2t. We also note that the hybridization term 
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Fig. 2. Examples of the energy band structure. The lower band becomes completely flat for e/ = V^jt - 2t. 

becomes zero for k - n, because of the form 2V cos{k/2) in the momentum space. Hence eigenvalue 
of Eq. (2.6) becomes -2t and e/. Figure 2 shows the flat band structure together with a case of sUghtly 
shifted €f for ? = 1 and V = 0.5. 

3. Variational Monte Carlo Method 

It has been proven in the case of equal number of electrons and unit cells that the ground state 
of this kind of models is ferromagnetic, provided the flat band condition is satisfied^^ or nearly so.^^ 
Ferromagnetism in two-band models near the flat band condition have already been discussed in great 
detail by a variety of methods. ^^'^^ In the preset work, we use a variational Monte Carlo (VMC) 
to investigate the ferromagnetism in CeRh3B2. In contrast with previous VMC work^^ for two-band 
ferromagnetism, however, our O-VMC takes account of spin dependence of effective hybridization 
and local level. As a result, our O-VMC has an advantage to provide the efi'ective band picture, which 
allows an intuitive explanation how the particular value of polarization is stabiUzed. This picture is 
especially useful for the case where the flat band condition is not satisfied strictly. Hence it is expected 
that the present work should bring new insight into the ferromagnetism near the flat band condition. 

We use an optimization technique^^ in the variational Monte Carlo (O-VMC) method. This tech- 
nique allows us to introduce five variational parameters. To consider a ferromagnetic state, we in- 
troduce spin dependent effective hybridization parameter Vg- and effective 4/ level e/o- Using these 
parameters, we construct the variational state as follows: We start from Eq. (2.6) with U - .By 
replacing V and e/ by VV and e/o- respectively, we obtain 

(ij)(T io- 

ia [ij]cr 

We diagonalize "Ho and derive the upper and lower hybridized bands which depend on the variational 
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parameters as 

Eka± ^\{^^k + ~efa ± V(^;fe-e/<r)2 + 16V^cosH/2j , (3.2) 
where ek = Itcosk. The creation operator of the upper (lower) hybridized band state aj^ i^la-^ is 



written as 



«L = ''kadl^ + vk^fl, (3.3) 

= -"kcrdl + (3-4) 

where uj^ and v^^ are weights of the 4 / and d states in the lower band: 

~ Elccr- ,„ 
Mfar = , (3.5) 
V4V2cOs2fc/2 + (6;t-£^^_)2 

2Vcosfc/2 

^ , ^ , =:• (3-6) 

V4y2cos2fc/2 + (6;t-£;t^_)2 

For a ferromagnetic state, we take the number A^^ of up spins larger than the number 'Ni of down 
spins without loss of generality. The total electron number 'Ng is given by A^g = /V^ + //j^. We first 
prepare an noninteracting state to construct a variational wave function. With A^-p < L, the magnetic 
state is given by 

where |0) is the vacant state, and No- over the product symbol means the number of /:'s involved. On 
the other hand, with A^^- > L for both spins, electrons occupy not only the lower band but also the 
upper band. Hence the noninteracting magnetic state is given by 

a- k k 

We have also considered a more general case where the upper hybridized band of up spins is partially 
occupied even though the lower hybridized band of down spins is not full. Starting from one of these 
noninteracting states, we construct variational wave functions by operating the Gutzwiller projection: 

^-W- H"fP' (3.9) 

where is a variational parameter to restrict double occupancy. \ir] = 1, any state with double oc- 
cupation is excluded, which is consistent with U = oo. On the other hand, if fj = 0, the Gutzwiller 
projection operator becomes the identity operator P = I, which corresponds to the case of U = 0. The 
variational wave function |*F) with the Gutzwiller projection is given by 

m^pm, (3.10) 

where the variational parameters are included in both |0) and P. 
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Fig. 3. Energy as a function of magnetization in the case of L = 120. 



4. Numerical Results for Ferromagnetism 

We define the magnetization m per site asm - (N-^ - N^)IL, and the density as Ug - (N-^ + N^)IL, 
where L is the total number of unit cells. We take V - 0.5 and either ef - -1.8 or e/ = -2.2 to 
determine the phase diagram in the plane of Ug and U. Figure 3 shows the energy per unit cell as a 
function of magnetization with Ug - \ and e/ = -1.8. We have taken the system with size L - 120 and 
the anti-periodic boundary condition. The minimum of the energy is located at m = Ug. In this case, 
the ground state is a ferromagnetic insulator and is fully polarized. We have done similar calculation 
taking L - 40 and found almost identical result with that in Fig.3. Therefore in the following we 
present results with L = 40 for various choices of Ug and JJ. 

We classify the ferromagnetic ground state into three kinds: 

(i) The fully polarized state where all spins align. This state appears ding - \. 

(ii) A band ferrimagnetic state which appears with 1 < < 2. Here the lowest effective band is fully 
polarized, and the next band is polarized in the opposite direction. Hence we have m-\-{ng-\)- 
2 - Hg. This state is also called the complete ferromagnetism in the literature.^^ 

(iii) A ferromagnetic state where the lowest band is fully polarized and the effective higher bands have 
no polarization. The magnetization of this state is m = 1, and appears in the region with high JJ and 
low electron density of the phase diagram for e/ = -2.2. 

(iv) Other ferromagnetic states which generally have m > 2 - n^. In this case the lowest band is fully 
polarized, but other bands are partially filled. 

Figure 4 shows phase diagram with L = 40 for (a) e/ = -1.8 and (b) e/ = -2.2. It turns out that 
the ground state can be ferromagnetic even in the region away from Ug - \. For ef - -2.2, the region 
with m = 1 is larger than that for = -1.8. Additionally the paramagnetic region expands for smaller 
U. 

To give a physical explanation of the ferromagnetism, we derive the effective band structure in the 
ferromagnetic state. First we consider the case of = 1.0 and U - oo, which is given in Fig. 5. The 
optimized variational parameters are given as - 0.5, e/| - -l.S, - 0.243 and e/^ - -1.52. 
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The lower up-spin band is not affected by the Coulomb interaction, while the lower down-spin band is 
uplifted. As a result all electrons occupy the lower up-spin band. 

Next we derive the effective band for = 1.3 and U = 0.5, which is given in Fig. 6. The 
lower down spin-band is affected by the Coulomb interaction and is uplifted as inne = 1 .0. Since the 
electron number is larger than the full occupancy of the lowest band, the extra electrons occupy the 
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Fig. 4. Ground state phase diagram in the plane of and U for (a) e/ = -1.8, and (b) e/ = -2.2. The radius 
of each circle represents the magnitude of magnetization. The dashed line in (a) separates the regions (ii) 
and (iii) as explained in the text. 




Fig. 5. The effective band for rig - l.O and U - oa. The circles represent occupied states of up spin electrons. 
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Fig. 6. The effective band for ng -1.3 and U - 0.5. The circles (squares) represent the occupied states of up 
(down) spin electrons. 

lower down-spin band. Thus the total magnetic moment has the magnitude m = 2 - rig. Namely the 
band ferrimagnetic state is realized. Since the parameters are close to the flat band condition, one may 
interpret this ferromagnetism as a generalized flat band ferromagnetism. 

On the boundary between (ii) and (iii) of the ground state phase diagram, the magnetization 
changes discontinuously. To clarify the behavior near the boundary, we compute the energy as a func- 
tion of the magnetization for U =1.5, 2.0, and 3.0. This result is shown in Fig. 7 . For U = 1.5, the 
energy has the minimum at m = 0.1 (Point A). As U becomes larger, another local minimum appears 
in the vicinity of m = 0.8 (Point B) The point B becomes the absolute minimum for U > 2.0. Hence 
the magnetization changes discontinuously at U ~ 2.0. 

We also derive the eff'ective band structure at point A and B for ?7 = 2.0. For U = 2.0, point A 
and B have almost the same energy. These results are shown in Fig. 8. At point A, the ground state 
is ferromagnetic which is connected to the flat band ferromagnetism. The effective lower band shifts 
to the higher energy compared to the one for U = O.ln the center of the lower energy band, the 4/ 
component is dominant and there appears almost localized states. These states become unstable as U 
becomes larger because the double occupation is not excluded completely. 

On the other hand, at point B, the electrons occupy not only the lower hybridized bands but also 
the upper band of up spins. The effective bands which consist mainly of 4/ orbital splits into up and 
down spin bands so as to reduce the double occupation. Thus the discontinuity of the magnetization 
is caused by the transition between the flat band ferromagnetism at point A and the almost local- 
ized ferromagnetism at point B. However, in reality some antiferromagnetism may also appear in the 
large U region. Since our variational wave function does not involve the antiferromagnetic correlation 
adequately, we cannot determine the most stable magnetic state by our approach. 

We have also derived the effective band for ey = -2.2, = 1.3 and U = 4, which is shown in 
Fig. 9. The Fermi level is located in the effective upper band. We find the optimum state such that the 
lower hybridized band of down spins and the upper hybridized up spins have the same Fermi wave 
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Fig. 7. Energy as a function of magnetization for U =1.5, 2, and 3 with L - 40. 




Fig. 8. The effective bands for rig - 1.9 and U - 2 with (a) w = 0.1 corresponding to the point A in Fig. 8, 
and (b) m = 0.8 corresponding to the point B. 



number. Namely, there is no polarization at the Fermi surface, and the value of kf agrees with those 
of the localized 4/ electron model. In contrast with the localized model, however, the fully polarized 
band involves the 4d component. Since we have taken the system size L = 40 in the calculation, we 
cannot exclude the possibility that there is a very weak polarization at the Fermi surface. 
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Fig. 9. The effective band for e/ = -2.2, tie - 1.3 and U - 4. The circles (squares) represent the occupied 
states of up (down) spin electrons. 



5. Application to CeRh3B2 

We compare the ferromagnetism of this model with experimental results for CeRh3B2. For this 
purpose we use the following parameters: f=0.34eV, ?7=7eV, ef = -0.714eV, V=0.24eV and rie = 1.1. 
Here t is determined by comparison with the band structure calculation.^ It is difficult to estimate 
the hybridization between the molecular orbital and the 4/ orbital. On the other hand, hybridization 
averaged over the band is estimated by X-ray absorption spectroscopy (XAS) to be between 0.23 and 
0.4 eV.'^"^' We tentatively adopt the value 0.24eV for V. The value ne = 1.1 is determined by the 
volume of the Fermi surface which is derived by the band structure calculation.^ Then we find that the 
ground state shows the complete polarization in the lowest band, and almost no polarization in other 
bands. The important features are as follows:: 

(i) Anisotropy of the magnetization 

The strong anisotropy of the magnetic moment can be explained by taking the crystal field state as 
= +1/2. In reality, other components of may be contributing to the ground 4/' level, and may 
influence the anisotropy and the magnitude of magnetization. 

(ii) Polarization of 4d electrons 

In the region close to the flat band condition, the lower band polarizes completely. Therefore 4d elec- 
trons also polarize significantly with strong hybridization. Because of the spin-orbit interaction of 4/ 
electrons, the magnetic moment of 4d electrons is anti-parallel to that of 4/ electrons. The magnetic 
moment per 4d electron is derived as O.I5/iB by our calculation, which is close to the value 0.18//B 
estimated by neutron scattering experiment. 

(iii) Reduced moment at Ce sites 

The total magnetic moment is estimated as 0.94/ib, which is larger than the experimental result 0.45;Ub. 
A possible source of difference is that the crystal field state assumed in the present work does not 
represent the actual electronic state. In order to reproduce the observed value of the moment, the 
expectation value of should be about 0.6 instead of 3/2 in the present model. 
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(iv) Fermi surface 

The effective band is shown in Fig. 10 where the Fermi level is located near the bottom of the effective 
upper band with the Fermi wave numbers kfa-ln = ±0.975. This situation corresponds to the case 
shown in Fig. 9. There is no polarization at the Fermi surface, and the value of kf agrees with those of 
the localized 4/ electron model. 

(v) Curie temperature 

Although our approach cannot discuss the Curie temperature directly, we can estimate the Curie tem- 
perature by the difference between the paramagnetic state (m = 0) and the ferromagnetic state. With 
the choice of the optimum magnetization m = I, the difference is calculated to be 826 K per site with 
the parameters used. This is surely the overestimate since our variational wave function in the param- 
agnetic state does not take proper account of correlation effects such as antiferromagnetic fluctuations. 
However, our result does show that the present mechanism of the ferromagnetism is consistent with 
the high Curie temperature of CeRh3B2. 

6. Conclusion 

In present work, we have applied the 0-VMC to study the mechanism of ferromagnetism near the 
fiat band condition. In the region both at and away from rig = I, we have found that the band ferromag- 
netism is stabilized. We have compared the ferromagnetism of this model with the experimental results 
for CeRh3B2. We have provided reasonable explanation for the anisotropic magnetization, reverse po- 
larization of conduction electrons, and the size of the Fermi surface probed by the dHvA effect. On the 
other hand, the total moment obtained by our model is much larger than the experimentally reported 
value. A possible source of difference is ascribed to our CEF state which may not be justified quanti- 
tatively. It is hoped that improvement of the present model will provide more quantitative explanation 
of the ferromagnetism in CeRh3B2. 
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